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The Navier-Stokes transport coefficients for a granular gas of smooth inelastic hard disks or 
spheres are determined from the inelastic Boltzmann equation by means of Grad's moment method. 
(NT The shear viscosity n, the thermal conductivity n and the new transport coefficient /i (not present 

for elastic collisions) are explicitly obtained as nonlinear functions of the (constant) coefficient 
of restitution a. The expressions of 77, k, and fi agree with those previously obtained from the 
Chapman-Enskog method by using the first Sonine approximation. A comparison with previous 
^ , results derived from Grad's moment method for two and three dimensions is also carried out. 
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I. INTRODUCTION 



The simplest model for a granular gas is a system constituted by smooth hard spheres or disks with inelastic 
collisions. The loss of energy in each binary collision is accounted for a constant coefficient of normal restitution 
a < 1. The essential difference from that for ordinary gases (a — 1) is the absence of energy conservation, yielding 
both obvious and subtle modifications of the usual Navier-Stokes (NS) equations for states with small spatial gradients 
of the hydrodynamic fields. For a low-density gas, the Boltzmann equation (BE) conveniently adapted to dissipative 
dynamics has been used in recent years as the starting point to derive the NS hydrodynamic equations. Therefore, 
\ assuming the existence of a normal solution for sufficiently long space and time scales, the Chapman-Enskog method 
[1] has been applied to determine the explicit forms of the relevant transport coefficients for dilute [2] and moderately 
dense [3] gases. As in the case of elastic collisions [1], the NS transport coefficients are given in terms of the solutions 
of a set of coupled linear integral equations which can be approximately solved by considering the lowest terms in 
a Sonine polynomial expansion. In this context, granular hydrodynamics follows similar steps as those made for 
ordinary gases. 

An alternative way of solving the BE is by means of Grad's moment method [4]. The idea behind Grad's method 
£j \ is to expand the velocity distribution function / in a complete set of orthogonal polynomials, the coefficients being 
the corresponding velocity moments of /. Then, the expansion is truncated at a given order k so that, the exact 
distribution function / is replaced by its truncated expansion f^ k \ When the approximation is substituted into 
the hierarchy of the moment equations up to order k, one obtains a closed set of coupled equations for the retained 
1 velocity moments. 

fT) , The inelastic Enskog equation was solved years ago by Jenkins and Richman [5, 6] by applying Grad's thirteen- 
Q\ ■ moment for a dense gas of inelastic hard spheres. Although the application of Grad's method to the Enskog or 
Boltzmann equations is not restricted to nearly elastic particles (a ~ 1), the results derived by Jenkins and Richman 
■ [5, 6] neglect the cooling effects on the granular temperature T due to the cooling rate. Given that this assumption 
t-H ' can only be justified for nearly elastic systems, their expressions for the NS transport coefficients differ from those 
CNJ . obtained in Ref. [2] from the Chapman-Enskog method for arbitrary degree of inelasticity. 

The goal of this contribution is to apply Grad's method to get the NS transport coefficients of d-dimensional 
granular gases at low-density. In contrast to the Jenkins- Richman theory [5, 6], my calculations account for the 
time dependence of the temperature coming from the inelastic cooling. As a consequence, the results are not limited 
a priori to weak inelastic systems. In particular, when only linear terms in the spatial gradients are retained (NS 
^ ■ approximation), the constitutive relations for the pressure tensor Py and the heat flux q are 



P i:i = P S l3 - r) i),t , + djUi - -6, X • U ) , (1) 



q = -reVT - fxWn, (2) 

where p = nT is the hydrostatic pressure, U is the mean flow velocity, T is the granular temperature, and n is the 
number density. In addition, 77 is the shear viscosity coefficient, n is the thermal conductivity coefficient and fi is 



'Electronic address: vicentegOunex.es; URL: URL:http://www. unex.es/eweb/fisteor/vicente/ 



2 



a new transport coefficient not present for elastic collisions. In contrast to the Jenkins- Richman theory [5, 6], the 
present results show that the expressions of 77, k and u obtained from Grad's moment method are the same as those 
obtained [2] from the Chapman-Enskog expansion in the first Sonine approximation. 

II. BOLTZMANN KINETIC THEORY FOR GRANULAR GASES 

We consider a granular gas composed by hard disks (d = 2) or spheres (d = 3) of mass m and diameter a. In the 
simplest model, the spheres are completely smooth so that the inelasticity of collisions is characterized by a (constant) 
coefficient of normal restitution a < 1. We also assume that the density n is sufficiently low so that the one-particle 
distribution function /(r,v,i) of grains verifies the (inelastic) BE. In the absence of external forces, the BE reads 

ft/ + v-V/ = J[v\f,f], (3) 
where the Boltzmann collision operator J[f, f] is 

J[vi|/,/] =o d ~ X \ dv 2 / d^e(^-g 12 )(^-g 12 )[ a - 2 /(r,vi;t).f(r,v^;t)-/(r,v i; t)/(r,v 2 ;t)] . (4) 

Here, d is the dimensionality of the system, <r is a unit vector along the line of centers, is the Heaviside step 
function, and gi 2 = vi — v 2 is the relative velocity. The primes on the velocities in Eq. (3) denote the initial values 
{v^Vj} that lead to {vi,v 2 } following a binary collision: 

vi = vi - i (l + a -1 ) (tr • gi 2 )CT, v 2 = v 2 + i (l + a _1 ) (er • gi 2 )S\ (5) 



The exact macroscopic balance equations for n(r, t), U(r, t) and T(r, t) follow directly from the BE (3) by multi- 
ying with 1, mv, and |mt) 2 and integrating over v. They are given by [2] 

D t n + nV-V = 0, (6) 



pDtUi + djPij =0, (7) 

D t T+^-(d iqi + Pij d Ui ) = -(T, (8) 
where D t = d t + U • V is the material derivative, p = ran is the mass density, 

Pij = J dv mViVjfiv) (9) 

is the pressure tensor, 

q = y dv |^ 2 V.f(v), (10) 

is the heat flux and 

is the cooling rate characterizing the rate of energy dissipated due to collisions. In the above equations, V = v U 
is the peculiar velocity. 



III. GRAD'S MOMENT METHOD 



As mentioned in the Introduction, Grad's moment method is based on the expansion of the velocity distribution 
function in a complete set of orthogonal polynomials (generalized Hermite polynomials), the coefficients being the 
corresponding velocity moments. In addition, to solve the (infinite) hierarchy of moment equations, the expansion 
is truncated after a certain order k and so, the above hierarchy becomes a closed set of coupled equations. In the 
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standard Grad's moment method, the retained moments are the hydrodynamic fields (n, U, and T) plus the irreversible 
momentum and heat fluxes (P^ — p5ij and q). In the three-dimensional case (d = 3), there are 13 moments involved 
in the form of the velocity distribution function /; hence this method is referred to as the 13- moment method. In the 
case of a general dimensionality d the number of moments is d(d + 5)/2 + 1. 

On the other hand, since we are interested in comparing the present results with those obtained for granular gases 
[2] from the Chapman-Enskog method, I'll include the full contracted moment of fourth order 



d(d + 2) 



AnT 2 



J dvV\f(V) 



d(d + 2) 



(12) 



The inclusion of the scalar field c to the thirteen moments of mass density, velocity, pressure tensor and heat flux vector 
will allow us to make a close comparison with the previous forms derived for the cooling rate and the NS transport 
coefficients [2]. Thus, the explicit form of the non-equilibrium distribution function for the fourteen moments can be 
written as 



771 2 711 f 111 \ f ^ 

/(V. - J,.V i{ l + -(P ij - P 6 ij) V i V j + —^(— 



d + 2 
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where 



d + 2 m,V 2 d(d + 2) 



(13) 



(14) 



is the local equilibrium distribution function. The coefficients appearing in each one of the velocity polynomials in 
(13) have been chosen by requiring that the basic fields, the pressure tensor, the heat flux vector and the fourth 
moment of the trial function (13) to be the same as those for the exact velocity distribution function /. It must 
be remarked that the fourth moment c characterizes the deviations of / from its Gaussian form in the homogenous 
cooling state. Note that, for the sake of simplicity, I have not included all polynomials of fourth order in the trial 
function (13). The inclusion of these moments would modify for instance the form of the cooling rate. However, as 
mentioned before, I want to offer here a theory with the same degree of accuracy as the one reported [2] by using the 
Chapman-Enskog method and so, only non-Gaussian corrections to the homogeneous cooling state distribution will 
be included in the Grad's distribution function (13). 

The cooling rate £ can be evaluated when one makes use of the form (13) for / in the definition (11). Neglecting 
nonlinear terms in P^ — p5ij , q, and c, the result is 



_ 2 7r< rf - 1 )/ 2 

Q -d r(f) ( 
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(15) 



The goal now is to determine the first order contributions to the irreversible momentum and heat fluxes and the 
fourth velocity moment c. Let us evaluate each quantity separately. 

A. Pressure Tensor 

The pressure tensor is defined by Eq. (9). In order to evaluate it, I multiply both sides of Eq. (3) by mViVj and 
integrate over velocity. The result is 



d t Pij + Pij V • U + U k d k Pij + Pk 3 d k U l + PktdkUj + j^dk(qiSjk + qjSik + qk&ij) = -^{Pij - P$ij) - (p$ij, (16) 
where use has been made of the relations 

J dv mViVjVkf{V) -> ^q-rj^Afc + QjSik + Qkfaj), J dv mViVjJ[f, /] ->• -v v {Pij - P$ij) - (pSij, (17) 



where [2, 7] 



(18) 
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and 

8 TT^" 1 )/ 2 . . [Y 

v o = ~, — t: t-tt — ncr \ — (19) 

d + 2 r(|) Vm v 7 

is an effective collision frequency related to the NS shear viscosity rjo of an elastic gas (^o = p/Vo)- Note that 
nonlinear terms in c, — p<5jj and q have been neglected in the collision integral (17) when / is replaced by its 
Grad's approximation (13). 

We are interested in the solution of Eq. (16) in the NS approximation, namely, when the pressure tensor and the 
heat flux are given by the constitutive equations (1) and (2), respectively. Thus, in order to solve Eq. (16), we need 
to make use of the balance equations (6)-(8) up to first order: 

d t n -> -U • Vn - nV • U, d t U -+ -U • VU - p _1 Vp, d t T -+ -U • VT - -|-pV ■ U - (T. (20) 

an 

The shear viscosity r\ can be easily obtained from Eq. (16) when one takes into account Eqs. (1) and (20). The 
corresponding equation for rj is 

(dt + Vr,)V=P, (21) 

where the time derivative dtrj must be evaluated to zeroth-order in spatial gradients. From dimensional analysis 
rj oc T 1 / 2 and so, dt"r\ — \f]dt InT = — \C,r\- With this result, the solution to Eq. (21) can be written as 

'-j^V (22 » 

where 770 = p/vo, is* = v ri /vQ, and = (,/vq. The expression (22) for 77 agrees with the one derived [2] from the 
Chapman-Enskog method. 

B. Heat Flux 



The evaluation of the heat flux follows similar mathematical steps as those made before for the pressure tensor. In 

m 
2 



order to determine it, one multiplies both sides of Eq. (3) by ^V 2 Vi and integrates over velocity. The result is 



d t Qi + \pd t Vi + PijdtUj + ^ll^SijdjipT) + d 3 v - {^y'Ph - + • u + u • v 9i 

djU k {q l 5 lk + qjSik + q k Sij) + U • VUjPlj + q-VU t + ^pU ■ VU l = -u K q u (23) 



d + 2 J " J " J * J ' J J " ' 1 2 

where use has been made of (17) and the relations 



/ dv I V 2 W(V) -+ dJ ^ V ^n + - p i^T-Pa P*<i) ■ / dv|y 2 VJ[/, /] -+ -^q, (24) 



where 

1 + a 



d-1 3 . , , 4 + 5d-3(4-rf)a ' 
_ + _( d + 8)(l-a) + 777^^ 



m. (25) 



As before, I have only retained linear terms in c, Pjj — p8ij and q in the second identity of Eq. (24) . In the NS 
approximation, Eq. (23) becomes 

dtqi ^r- P diP + i£r ( l + i) di{pT) = ~ VKqi - (26) 

Dimensional analysis shows that k oc T 1 / 2 and u oc T 3 / 2 . Thus, according to the constitutive equation (2), the time 
derivative of the heat flux can be written as 

dm = 2Cnd t T + (C^ + l^j dm, (27) 
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where use has been made of the relation 



di{d t T) = -diiCT) = -(Tdi Inn - -QTd, InT. 



(28) 



Substitution of Eq. (27) into Eq. (26) allows one to identify the explicit forms of the transport coefficients n and fi. 
They are given by 



d-1 1+c 

K = ; Kq- 



d %*-2C*' 



(29) 



KqT 



» A * , d - 1 



where 



d(d + 2) r)o 
2(d- 1) m 



(30) 



(31) 



is the low density value of the thermal conductivity of an elastic gas, v* K = v K /vo, and n* = k/ko- As in the case 
of the shear viscosity, Eqs. (29) and (30) for k and /z, respectively, agree with those previously derived [2] from the 
Chapman-Enskog method in the first Sonine approximation. 



C. Fourth Velocity Moment 



To close the complete determination of the NS transport coefficients, it still remains to evaluate the contracted 
fourth moment c defined by Eq. (12). To get it, I multiply now both sides of Eq. (3) by V 4 and integrates over 
velocity. In the NS order, one gets 



d(d + 2) (l + |) m- 2 (d t (pT) + ^PTV • U + U • V(pT)) = J dvV 4 J[f, /]. 



(32) 



Equation (32) can be simplified when one takes into account the balance equations (20) and the relation (which only 
applies to linear order in c) [8] 



J dvV-J[/,/M-^g(l-« 2 ) 
In dimensionless form, Eq. (32) can be written as 



< i+ i + « 2 + f(i< 1M+39+lto2 > + S 



Vq. 



Ad\ 



(l + |) C = (1 - « 2 ) d + \ + + °- (A(10d + 39 + 10« 2 ) + | 



d-1 



(33) 



(34) 



The solution of Eq. (34) can be easily obtained when one uses the expression (15) for the cooling rate (. Retaining 
only linear terms in c, the final form of c is 



32(1 - a)(l - 2a 2 ) 



9 + 24d - a(41 - 8d) + 30(1 - a)a 2 ' 
Equation (35) coincides with the one derived by van Noije and Ernst [8] in the free cooling state. 



(35) 



IV. COMPARISON WITH PREVIOUS RESULTS 



The results derived in the previous Section from Grad's moment method have shown that the explicit forms of 
the NS transport coefficients are the same as those previously obtained in Rcf. [2] by using the Chapman-Enskog 
method. This shows the equivalence between both methods in the NS regime. On the other hand, as said in the 
Introduction, Grad's 13-moment method [4] was already used many years ago by Jenkins and Richman to determine 
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the NS transport coefficients of a dense gas of inelastic hard disks (d = 2) [5] and spheres (d = 3) [6]. In the low-density 
regime and for smooth spheres, their expressions are given by 

Vm = -jz — Q 8 w-, : — r (d = 2), T) JR = — — - (d = 3), (36) 

(7 — Sa)(l + a) (3 — a)(l + a) 



kjr = 7— — r 2 ^° 15 , -r (d = 2), K m = — — r-, 2 ^^ vT ( d = 3 )' WR = 0. (37) 

(1+a) (l-a)J (1 + a) [l + ff(l-a)J 

A more recent derivation [10] of the NS transport coefficients for hard spheres has been also carried out from Grad's 
fourteen moment theory. It differs from the Jenkins and Richman approach [5, 6] by the inclusion of the contracted 
fourth moment c. In particular, the coefficient /i is different from zero but very small since it is proportional to c. 

Comparison between Eqs. (22) and (36) for rj and Eqs. (29)-(30) and (37) for k and /i show that the results derived 
in the present paper differ from those obtained by Jenkins and Richman [5, 6]. The differences between both works 
are due essentially to the assumptions made in Refs. [5, 6] since the latter authors neglect the time dependence of 
temperature due to collisional cooling (which is equivalent to take £* = in our expressions) and the non-Gaussian 
corrections to the homogeneous cooling state (which is equivalent to take c = in our expressions) . While the latter 
simplification is in general not relevant (except for quite extreme small values of a), the former hypothesis turns out 
to be important for finite values of dissipation (i.e., beyond the quasielastic limit a ~ 1). In particular, the coefficient 
\x vanishes in the Jenkins- Richman theory [5, 6]. 

In order to illustrate the quantitative differences between both theories, Fig. 1 shows the a-dependence of the 
(reduced) transport coefficients if (a) = rj(a)/r]o, K*(a) = k(o)/ko, and fi*(a) = n/x(a)/T«o for hard disks (d = 2) 
and spheres (d — 3). It is quite apparent that while the dependence of the shear viscosity on the coefficient of 
restitution is relatively well captured by the Jenkins-Richman theory, there are significant differences for the heat 
flux transport coefficients. In particular, the present theory shows that n* increases with decreasing a while the 
opposite happens in the Jenkins-Richman theory. The discrepancies are much more important for the coefficient fi 
since my results show clearly that fi ^ in contrast to the prediction of the Jenkins-Richman theory (/ijR = 0). The 
comparison with the results derived in Rcf. [10] are not presented here since the latter results slightly differ from 
those obtained in the Jenkins-Richman theory for spheres [6]. It is important to remark that the Chapman-Enskog 
expressions for r/* , k* and fi* (which are the same as those obtained here) compare quite well with numerical results 
[9] obtained by solving the BE by means of the direct simulation Monte Carlo (DSMC) method. 

In summary, Grad's moment method has been used to determine the NS transport coefficients of a granular gas 
described by the inelastic BE. As in a previous work [10], the present method differs from the conventional 13-moment 
method by the inclusion of the full contracted fourth moment c, defined by Eq. (12). In contrast to previous attempts 
[5, 6], the present theory takes into account the time dependence of the granular temperature due to cooling effects as 
well as the non-Gaussian corrections to the homogeneous cooling reference state. The explicit expressions for the NS 
transport coefficients derived here agree with those obtained years ago from the Chapman-Enskog expansion in the 
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first Sonine approximation [2]. This agreement shows the equivalence between both approximate methods to solve 
the BE for granular gases in the NS regime. The inclusion of more velocity moments (for instance, all the fourth 
degree velocity moments) would change the final results since there would be likely additional contributions to the 
cooling rate, for instance. As a future work, I plan to extend the present theory to moderate dense gases [3] by using 
the Enskog kinetic equation. 
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